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Convex Sets in Linear Topological Spaces 


ISRAEL HALPERIN, F.R.S.¢ 


INTRODUCTION 


|’ this paper we shall prove some theorems concerning convex sets 
in linear topological spaces. 

In §2 we shall show that if A and A, are non-empty, disjoint convex 
sets in a real-linear topological space and if at least one of them has 
inner points then there is a hyperplane separating them. This was 
proved by M. Eidelheit [3] for the more special case that the space is 
linear and normed and both A and A, have inner points 

From this result we deduce that if A is a convex set with inner 
points then through every boundary point of A there passes a sup 
porting hyperplane; more generally if R is a linear set containing no 
inner points of A then there is a hyperplane // containing R such 
that A lies on one side of /7. This was proved for linear normed spaces 
by S. Mazur [4], the first half of this result having been obtained 
previously by G. Ascoli [1] for separable, linear, normed spaces 

Aside from greater generality, our proofs are more elementary than 
those of Ascoli, Mazur, and Eidelheit and do not use Minkowski 
functionals or the Hahn-Banach extension theorem. (The simplified 
proofs for Eidelheit’s theorem due to WKakutani [5] and Botts [2] 
require the theorem of Ascoli-Mazur 

In $3, we consider linear topological spaces with either the real or 
complex numbers as scalars. We repeat the argument of Ascoli and 
Mazur to obtain tke theorems (proved by them for the real-lineat 
normed case) that every closed convex set is weakly closed and that 
if x, ~ x weakly then x is in the closure of the set of all averages 
CiX1 +... + Cm¥m With m  S ., all c, > O and the sum of the 
cn equal to 1; however, for this result we require the stronger re 
striction on the space that every open set containing the origin should 


contain a convex set having the origin as inner point. 
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2. SEPARATING HYPERPLANES 

Let / be a linear space topologized by a family of open sets such 
that: the set-union of an arbitrary collection (possibly empty) of 
open sets 1s an open set, the set-intersection of two open sets is an 
open set, and each point of & is the intersection of all open sets con- 
taining it. We shall suppose (i) x + y is continuous in x for fixed y 
and (ii) cx is continuous in the pair c, x. It follows that whenever 0 
is in the closure of a set of x then O is also in the closure of the cx 
where ¢ may depend on x providing that A = sup/c| is finite; 
indeed, for arbitrary integer p each cx may be expressed as m2>?x 
+ rx with O<r, < 27”? and m one of the finite set of integers 
numerically less than 2?(2 + A), and hence by choosing p sufficiently 
large it can be seen that in every open set containing O there is at 
least one cx 

If S,, S» are subsets (possibly single points) of # then S$; + S» shall 
denote the set of all x, + x» with x; in S; and x» in So, and cS the set 
of all cx with x in S; J(S,, S.), the join of S; and S», shall denote the 
set of all cx; + (1 c)xo with x, in S;, x. in S., and 0 < c < 1. Our 
conditions on / imply that x + N and cN are open sets for fixed x 
and c # O respec tively if N is an open set. 

\ set A is called convex if A J(K, K). It is easily verified that 
T(S,, So) is convex if 5S), Se are both convex and that the set union of 
an increasing family of convex sets is again convex. 

Our conditions on / ensure that if AK is convex then its closure K 
is also convex; for if K C K, for some convex K, C K and if x is in 
K then J(x, Kj) is again convex and K C K, C J(x, Ay) C K and it 
follows that AK is contained in a maximal such A, which must then 
coincide with K. 

lor the rest of this section we shall suppose that the scalars are the 
real numbers. 


} 


I(x, y) shall denote the set of alltx + (1 — ty with — © <¢ « x 


I(x) shall mean /(O, x). A set /7 is called linear if 77 contains /(x, y) for 


all x, yin //; [7 is called a hyperplane if it is closed, linear, and different 
from £ and the only linear sets which contain /7 are // and 
Suppose // is a hyperplane. Choose any xo in 1, any yo not in I 
and set u Vo vo. Then the linear set /7 + /(w) must coincide 
with & and hence every x in / can be represented in the form h + tu 
with / in J/. Clearly 4 and ¢ are determined uniquely by x. Let 0 have 
the representation hy + tou and for every x = h + tu set f(x t to. 
hen f(x) is a real-valued linear functional on / and our conditions on 
FE. together with the closure of 7/7 ensure that f(x) is continuous. (If 
f(x) were not continuous we could deduce that for some ¢ > 0 and 
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some fy, t; and some set of A, ¢, | t t,| > € but A; + fu Is in the 
closure of h + tu; hence 0 is in the closure of the set of h h, + 
(t — t))u; then, since | t — 4, ! < € ', 0 is in the closure of the set ot 
(t — 4h) "(h — hy) + vo xo and finally yo is in the closure of /7/.) 
H consists precisely of the x for which f(x) = — fo. Each of the sets 
of x for which f(x) < — to, and f(x) > ty respec tively, is called 
a side of IH and also a half-space. It is easily verified that the two sides 
of /7 are uniquely determined by // and that if y is an inner point of a 
set S lying on one side of /7 then y is not in // 

Conversely, if /7 is a closed linear set with J] # E IT + I(u) tor 
some u then it is easily verified that // is a hyperplane. Also, if f(x) is 
any real-valued continuous linear functional not identically zero and 
cis areal number then the set of x for which f(x) = ¢ is a hyperplane 

We shall say that two sets S; and S» are separated by the hyperplane 
7 if they lie on different sides of /7, equivalently if, for some continuous 
real-valued linear functional f(x) and some c, f(x) cif and only if 4 
is in HJ, f(x) < c for all x in S; and f(x) > c for all x in S 

A hyperplane /7 is called a supporting hyperplane of aset S if S lies 
on one side of // and // contains at least one boundary point of S 


LEMMA. Jf Cis a closed set with inner points and its boundary IT 1s 
a non-empty linear set then Il is a kyperplane and C ts one of its sides. 
Proof. (i) By a translation, if necessary, we may suppose that the 

origin O is in JT 

(ii) Since C is closed and its boundary is a linear set it 
follows that C is convex. Hence if u is in C then ¢u is in C for all t > O 
and h + tu (1 4 (1+ 2)7'h + t(1 + t)~'%) ts in C for all h in 
HT and t > 0 

(ii) If «and —w are both in C then uw is in //. For otherwise 
there would be an open set N containing u and contained in C. It 
would follow that /(.V, —1) is contained in C and hence 4(.NV + (—u 
which is an open set containing the origin, is contained in C, contra 
dicting the fact that the origin is a boundary point of C 

(iv) If wo is a fixed point in C but not in // ther C consists 
precisely of all h + tuo with h in /7 and ¢ > 0. Indeed, it v in C were 
not of this form it would follow that for every t > 0, v — tuo is not 
in /7. Since v would be an inner point of C it would follow that for 
every t > 0, v — tuo is in C, t"'(v — tuo) t—'4 uy is in C. When 
t— © we obtain — wp is in C, assumed closed. By (iii) it would follow 
that wo is in /7 contradicting our assumption 


(v) With the uo of (iv), JZ + 119) I. since’it is a linear 


set having inner points. Since // is closed this implies the lemma 
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THEOREM |. Suppose that K and kK, are non-empty disjoint convex 
ets and that K ts open. Then there is a continuous linear f(x) and a « 
such that f(x) < c¢ for allx in K and f(x) > c for all x in K,, so that the 
hyperplane f(x) = c separates K and K,. 

Proof. (i) Since the union of an increasing family of convex sets, 
each disjoint from K, is also a convex set disjoint from K we may 
suppose that A, is not contained in any other convex set disjoint 
from A. Then A, will be closed since the closure of A, is convex and 
disjoint from the open set K. 


(ii) If wisnot in A, then J(u, Ky) # Ay. But the set J(u, A,) 
is convex and contains A, so it must intersect K. Thus cu + (1 
must be in A for some x, in A, and some 0 < ¢ < I. 


> 


C)X, 


(iii) If x, isin A, and uw is arbitrary then at least one of the 
sets: all x; + tu, t > 0 or all x, tu, t > 0 is contained in A,. Other- 
wise there would exist “,; and m#, neither in A, and with x; in J(1, ue) 
Then by (ii) there would be y,, 2; with z,; in Ay, y, in K and y,; contained 
in J(m, 2); similarly for ts, yo, 2. Then J(;, ye) would intersect the 
triangle of x), 2 and 2, that is, J(%,, J(%, 22)), in a point w which 
would be in both A and A,, a contradiction. 


(iv) If x, is in A, but uw is not in A, then for every t > 0, 
Kr, + U(x; wu) is an inner point of Ay. For there is an open set N 
containing u and disjoint from Ay. Then x; + ¢(x; — y) is in A, for all 
t > 0 and all y in N, and for fixed ¢ > 0, x1 + U(x, + (—1I)N) is an 
open set containing x) + (x, u). 


(v) Since A, is a convex set it now follows that the boundary 


of A, is a linear set and that A, is a closed set with inner points. The 


theorem then follows from the previous lemma. 


COROLLARY |. Jf Ko and Ky, are non-empty convex sets and K has 
inner points, then K and K, can be separated by a hyperplane whenever 
K, ts contained in the closure of a convex set K’ such that K’ contains no 
inner point of K, in particular if K, contains inner points but K and K, 
have no common inner points. 

Proof. This follows from the fact that if a convex set has innet 


points it is contained in the closure of the convex set of its inner points. 


COROLLARY 2. Jf x is a boundary point of a convex set K possessing 
inner points then there exists a supporting hyperplane of K which 
contains x 


Proof. This follows from Corollary 1 if x is taken as Ay 


PHeOREM 2. Suppose that K is a convex set with inner points and 
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that R is a non-empty linear set containing no inner points of K. Then 
R 1s contained in a hyperplane H such that K lies on one side of H. 

Proof. Replacing AK by the convex set of its inner points, we may 
assume that A is an open set and that A and R are disjoint. 

Set J = J(R, K). Then J is a convex set with inner points. More- 
over no inner point of J is in R. Indeed, if x in R were an inner point 
of J we could choose any x» in AK and deduce that for some ¢ > 0, 
x + t(x — xo) = 1 isin J but not in R and x is in J(xo, 1). Then 
would be in J(x, y) for some x; in K and y in R. This would imply 
that /(x, y), all in R, intersects J(xo, x;) all in A, contradicting the 
assumption that A and R are disjoint. 

Corollary | to Theorem | now implies that J and R can be separated 
by a hyperplane //. Since R is contained in J, R must lie on both sides 
of HT and hence is contained in J7. 


THEOREM 3. Jf FE is such that every open set containing 0 contains 
a convex subset with 0 as an inner point, then if K is a non-empty closed 
convex set and y ts not in K there is an open set N containing y such that 
K and N can be separated by a hyperplane H. (Necessarily, y is not in H.) 

Proof. By our assumptions on £ there will be an open set N con- 
taining y with N contained in a convex set A, such that AK and A, are 
disjoint. Since A, has inner points it follows from Corollary 1 to 
Theorem 1 that AK and Ky, a fortiori K and N, can be separated by a 
hyperplane. 


COROLLARY. For E as in Theorem 3, every closed convex set different 
from EF 1s the intersection of the half-spaces which contain 1t. 
3. WEAK CONVERGENCE 


If & is a linear space with a topology as specified in §2 we shall say 


that a sequence x, in / converges weakly to x if f(x,) converges to f(x 


for every continuous linear f(x) on / which is real or complex valued 


according as the scalars for / are the real or complex numbers. A set Z 
contained in £ is called weakly closed if x isin Z whenever a sequence x, 
in Z converges weakly to x. 

If the scalars are the complex numbers then / may be regarded also 
as a linear space /’ with the real numbers as scalars and with the same 
open sets as £. If f(x) is complex-valued and linear and continuous on 
E, then clearly (x) h f(x) is real-valued, linear and continuous 
on ££’ and conversely for every real-valued, linear, continuous (4 
on EF’, f(x o(x i@(1x) is complex-valued, linear and continuous 
on F:. 
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THEOREM 4. Jf every open set containing the origin has a convex 
subset with the origin as inner point, then: 

(i) Every closed convex set 1s weakly closed. 

(ii) If x, converges weakly to x then x 1s in the closure of the set of all 
averages (\X; +... + CmXm with m = 1, 2,...,allcen > Oand the sum 
of the c, equal to 1. In particular, if E. satisfies the first countability axiom 
of Hausdorff there will be a sequence of such averages which converge to x 
in the sense of the open-set topology. 

Proof. We may suppose that the scalars are the real numbers since 
¥, converges weakly to x for complex scalars if and only if this is true 


for real scalars 


If now K is a closed convex set and x, in AK converges weakly to y 


and y is not in A then Theorem 3 implies that there is a hyperplane // 
separating A and y with y not in //. Then there is a continuous real- 
valued linear f(x) and a real number c such that f(x) < c¢ for all x in 
K but f(y) > c. This contradicts the fact that f(x 

to f(y). 


‘n) Must converge 


(After this note was accepted for publication the writer learned 
that an elementary proof of the Mazur-Ascoli theorem had been found 
earlier by M. H. Stone and others, see V. L. Klee, Jr., Convex sets in 
linear spaces, Duke Math. J., vol. 18 (1951), pages 443-66, 875-83 
for references. ) 
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Bessel Expansions of the Confluent Hypergeometric 
Functions* 


r. E. HULL, C. A. SWANSON, and D. A. TRUMPLER 


Presented by R. D. JAMES, F.R.S.C. 


INTRODUCTION 


A’ initial or boundary value problem associated with a dif 
ferential equation can be reformulated in terms of an integral 
equation in a variety of ways. When the problem cannot be solved 
exactly, such an integral equation is usually the starting point of any 
successive approximation procedure used to generate the solution. 
It is also natural to begin with such a reformulation when investi- 
gating the asymptotic behaviour of the solution, or, in the case of 
boundary value problems, when investigating the convergence of an 
eigenfunction expansion or the asymptotic behaviour of eigenvalues 
and eigenfunctions. 

We shall first discuss briefly the reformulation of a general initial 
value problem in terms of an integral equation. We then adapt the 
method of successive substitutions to a special case and obtain the 


multiplication theorems for the Bessel functions. An extension of 


the idea leads to expansions of the confluent hypergeometric functions 


in series of Bessel functions. 

We shall include references to the other applications mentioned 
above and to other expansions of the confluent hypergeometri 
functions which have been found by different methods 


THE INTEGRAL EQUATION 


Suppose that after dividing through by p(z), the coefficients in the 
differential equation 


d dy 


*This investigation was carried out while two of the authors (C.A.S. and D.A.-1 
held National Research Council Scholarships 
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are analytic in the finite z-plane except perhaps for poles at the origin. 
Then, for « # 0, 


(2) y(z) = n(z) + bue(z) + ¢ | {vi(x) ve(z) — v4 (2) vo(x)} 
e € 
r(x) v(x) dx, 


where 7, ve are independent solutions of 


d dy 
(3) ds ( p12) 2) + q(z)1 0 


and where a, 6 are arbitrary constants while 


’ 


c = 1/p(x) {u(x) v2'(x) V(x) v,'(x)} 


The transformation from (1) to (2) links the solutions of (1) with 
those of the associated equation (3). 


The transformation can be generalized to apply to n-th order 
equations. For example, with suitable restrictions on the coefficients, 
the differential equation 


n—1 


n—1 
(n) i) (i) 
(1’) y(z) + De ails) v'P(z) = DY rilz) y“(z) 
jam () j= 

can be transformed into 


n 


(2’) y(z) = >, a, 2;(z) }- | K(z, x) y(x) dx, 
im] Ve 
where the v,; are independent solutions of 


n—I1 
‘ (n = (a 
(3°) v” (2) + Dd giz) v'" (2) = 0. 
K(z, x) can be found by the method of variation of parameters and it 
turns out to be 


K(z, x) 


l 
i 


Bis 
hi ) V(z, x) r(x) ¢ 


where 
v(x) 


Vv, (x) 
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The transformation from (1) to (2) was first used by Liouville [10] in 
the special case when p(z) 1, g(z) = p*; the resulting integral 
equation was then used to study the asymptotic behaviour, for large p, 
of the eigenvalues and eigenfunctions determined by (1) after the 
imposition of certain boundary conditions. More general problems 
of the same nature have been considered by Birkhoff [2; 3] . 

Many papers have been devoted to a study of the asymptotic 
behaviour of solutions of differential equations. A majority of these 
papers have been based on a reformulation of the problem in terms 
of an integral equation and the pertinent theory has been developed 
chiefly during the past twenty years. Recent developments are given, 
for example, by Langer [14] and Cherry [7]. A transformation of an 
n-th order differential equation has been considered in detail by 
Miller [15]. 

It should be noted that the terms on the right side of (1) or (1’) 
need not be linear in y and its derivatives. If these terms are non- 


or 


linear, but relatively small, the form (2) or (2’) is useful in finding 
successive approximations to the solution if the problem is an initial 
value problem. A trivial but interesting example of this fact arises 
in relativity theory where the path of light past the sun depends on 


the solution of [8] 


du 


n bs 
dé + u 3mu',u(O0) = R (0) = 0 


where m is small. (Other methods are more direct when only the 


so-called limit cycle solutions of a non-linear problem are required 


MI ULTIPLICATION THEOREMS FOR BESSEL FUNCTIONS 


In the above mentioned studies of asymptotic expansions, the 
usual procedure has been to choose the integral equation (2) or (2’) 
so that the integral contained in it is of small order compared to the 
remaining terms, as some parameter approaches infinity. The re- 
maining terms on the right side are then asymptotic to the desired 
solution y(z). Apparently Ikeda [11] was the first to use the method 
of successive substitutions in (2) to obtain complete and convergent 
expansions of the solution in special cases. His results are not entirely 
correct; we shall therefore now treat his most important example in 
order to illustrate the procedure 

The differential equation for J,,(az) (vy not an integer) can be 
written in the form (1) as 


d (. dy 
(4) dz ~ ag 
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In this case », = J,(z) and v J_,(z), and (2) becomes 
(5) y(z) =a J,(2) + b J,(z) 


gli-—a) f°: 
a | {J_,(x) J_,(z) — J,(z) J_,(x)} x y(x) dx. 
Z2sinvmr oJ, ; 

Here and in later sections we need, for i = 0 


results* 


us 


(6a) | BE C2 ae ae Cy a e 
Z2sinvard, 


where 


(6b) A, 


$(2 + 1) sinvr 


Te 


(1 1) sin vr 


and 


a T 
(/a) > 
~ SILVWs 


where 


- , Tre 
(7b) ¢ . Fhe 
t+ l)sinmvarte 


Tre 


1) sinvr 


tf 
{ 
1(€) 1. |. 
, de 


*To avoid ambiguity in these results we will assume the s-plane to be cut, say 


along the negative real axis. The expansions obtained later will nevertheless hold 


without any restriction on arg 
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These results can be established as follows: it is a straight-forward 
matter to show that the left side of (6a) satisfies 


d du y ss 
, ‘ae so J n44(8). 
dz (. is) 1 (. 4 u a Junetd) 


The first term on the right of (6a) is a particular integral of this 
equation and the other terms are the complementary function. A, 
and B; are found by putting z = € in (Ga) and in the equation obtained 
by differentiating (6a). The proof of (7) is analogous. 

kor each a and 6, the method of successive substitutions applied 
to (5) will generate a series which converges to the solution of (5) for 
any finite z, e # 0. Assuming for the moment that we can rearrange 
the various contributions to the terms in this series, it is clear from 
(6) and (7) that we can write the solution in the form 


(8) y(z) =. Piz J v4.:(2) + > giz'J »~ 4S), 
0 0 
where p; and q; are constant coefficients. 
Proceeding formally, we now substitute (8) into (5); using (6) and 


(7) and equating coefficients of like terms, we then find that (5 
satished if p; and q; satisfy the recurrence relations 


is 


| P 
(9) : — _ 0,1,2 


2(i + 1) 2 l Byeres 
The values of po and qo depend on a and 6 and so ultimately on the 
particular solution of (4) which is being considered 

kor the solution to be J,(az), which is finite at the origin when 
R(v) > 0, we must take go = 0 in (8). Equating coefficients of 2” on 
both sides of (8), we obtain the value of po. From (9) the solution 


then formally reduces to 
(10) J (az) 


which we shall now show holds for all R(v) > O 
Since [20, p. 44] 


(2/2) 
(11) J,(%)| « exp 


(vy 7 1) tivo + 


where | vo + | | is the smallest of the numbers | v 
the series (10) is dominated, when R(v) > 0, by 


‘ » | 2 « 2 
(az/2) Z 
exp 


(vy + 1) 
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which obviously converges. The convergence is uniform with respect 
to z for z in any finite region, so our formal operation of integrating 
term by term with the aid of (6) is justified. It can be verified that 
the subsequent rearranging depends only on the absolute convergence 
of both (10) and the series obtained from (10) by differentiating. 
Requiring the coefficients of J,(z) and J_,(z) to vanish serves only to 
define a and b, and all the remaining coefficients vanish automatically. 
The convergence of (10) is also uniform with respect to v, and hence 
we can let v assume integral values. The formula (10) is therefore 
valid for all R(v) > 0. 

lor the solution of (5) to be J_,(az) (R(v) > Oand v not an integer), 


we find on equating coefficients of 2” on both sides of (8) that po = 0. 


I'(v) sin vr (:) ¢ 
1 i 


Using this result and (9) and putting y = J_,(az), we find, on equating 


further, near z (0), 


coethcients of (2/2) in (8S), that 


(vy) sinvr _, gow (—1) 01 a)" < ; 
a : ® : (vy + 2) sin (vy + 4) 2, 
w T= 


Tv 1. 


so that finally g a’ provided | | ai <-d. 


Using (11), we find that the series for J_,(az) is dominated by 


(2a) T'(v) sin vr Z —~il—a’l \Twt+z) 
: exp 2. 
WZ bivgo + 1! T% i! IT (y) 


which converges 1! 


\s before we therefore obtain, for R(v) > O and vy not an integer 
and | | 


Moreover, 


J (az) cos vr 


Sin vqT 


' Jya1(2) CoS ve ( 


Sin vw 


| 
‘) J y4i(z) cos (vy 4+ 
¥ sin (v 
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that is, 


(12) Y,(az) =a’), 


Equation (12) therefore holds, at least when R(v) > 0 and » is not 
an integer, so long as | 1 — a?| < 1 and g is finite and different from 
zero. That (12) holds under the same restrictions on @ and z even 
when » is an integer n, follows from the fact that 


V ni; (2) O((n + 14) !) 


(so that the series converges when v = ») and from the continuity of 
VY, as a function of v. 


Equations (10) and (12) are the so-called multiplication theorems 
for Bessel functions [20, p. 142]. 


I.XPANSIONS OF THE CONFLUENT HYPERGEOMETRIC FUNCTIONS 

We shall begin with the differential equation for the confluent 
hypergeometric functions My,,m(y) given by Whittaker [21, p. 337] 
dW ; a m*\ 


-.4 —(W =0, 


13 
; dy ' 1 : y 


where 2m is assumed to be net an integer. The change of variables 
and parameters 


(14) tky e. Wy) 2F(z), 2m 


transforms (13) into the form of (1): 


ee =) ( se oe 
dz \"ds/ * ; 16k 


where v is not an integer, so that 2, J (2), 0 J_,(z), and the 


v\elsys . ¥ 


corresponding integral equation for F becomes 


(15) F(z) a J,(z) + b J_,(z) 


1 a a i 
‘ Jy v) J y\Z) ‘ S). pit )t xX F(x) dx. 
32 k° sin a | 


Rearranging the terms in the series generated by the successive substi 
tution procedure can again lead to a solution of the form (8). On 
putting (8) into the integral equation (15) we find, however, that 
before each integration can be performed with the aid of (6) and (7), 


it is necessary to use a recurrence relation on each term of the series 
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Precisely, we rewrite (8) as 
: os yr 
F(z) = >> pa'| 2 

jax 0 z 

j a 

} qi 


im 


before substitution into the integral of (15). The relations (6) and (7) 
can then be applied. Proceeding formally, and putting Pp, Popi, 


q, Jodi, Po qo | we obtain, for 7 | Ae the relations 


(16 


i+3 «ps 2 
(4 32k 
\rguing as we did to establish the result (10) for J,(a@z), and returning 


to the original variables through (| 1), we obtain 


(18) Miam(y) Ake” "DP (2m 4 1)>. DA2V ky) Tony (2 ky) 


iO 


which is valid for all k, finite y, and for R(m) > 0; the p; are given 
by (16). Using Langer’s method, Taylor [18, Eq. (83)] has shown 
that the first term of this series is asymptotic to My,m(y) as| Rk — @. 

Phe corresponding result for A/;,,-n(y) with R(m) > O turns out 


to be an asymptotic expansion valid for large | Rk |. It is 


(19) M,.-m(y) ™ qo Gi(2V ky) 'T 2m (27 ky) 


for y finite and different from zero; the q, are given by (17). 
We shall first indicate how to prove that (19) is an asymptotic ex- 
pansion vaiid for large |k > , and then find go (see equation (21)). 


Substituting the expression 
20) F(z) : pod, pa'In, (z) 4 god, Giz’ J »-i(2) + f(z) 
t=Q iw 


into equation (15) and carrying out the integrations over the partial 
sums, we obtain an integral equation for f(z). Using the recurrence 
relations (16) and (17) and choosing a, 6 or po, Go so that the co- 
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efficients of J, and J_, cancel out, we finally reduce the equation 
for f(z) to the form 
) 


f(2) = 5642 [Pn-2 F Pui 


qo v i { . . 
: eee ae : tJ (e)J_»(s) J,(s)J_,»(x)} 
° L6R° lin: 32k sin var J. ee ys 


x f(x) dx 


Now, if J denotes the maximum of | f(x) | along the path of inte 
gration, then 


Ucg(|poi| At+iq| B+ MC)/|k 


where A, B, C are upper bounds for corresponding expressions in the 


preceding equation. Therefore 
M 


if | k | is chosen large enough that the denominator is positive. From 
(16) and (17), both Pi, Ji are Of k 1/3) so that 


lin+l 


M = po O(|kl te”) + go ORI 


moreover, the last term in each of the partial sums on the right of (20) 
is O( | k |~?"*). The series obtained from (20) is therefore asymptotic 
to the solution of (15) for large | R 

Equating coefficients of z’ and z~” in the F(z) corresponding to 


Mx,~m(y) [21, p. 338] and in the right side of (20), we obtain 


po = 0 


v1) /2 . . 
(4k) ~ go(sin vr r) >, (—1)'G.0'( + 2). 


Hence, returning to the original variables, we obtain an asymptotic 
expansion for the coefficient go in (19 


(21) 1/qo~ (4k) "8 (sin 2ma/mw) >> (—1) G0 (2m + i). 
Appropriate changes of variables in (18) and (21) lead to ex- 
pansions of, for example, the Laguerre and the Hermite functions. 


CONCLUDING REMARKS 


We have shown that the usual method of studying the asymptoti 
solutions of differential equations can, with some modification, be 


used to obtain complete expansions ot spec ial functions in certain 
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important cases. In a future paper we propose to adapt the method 


to the solution of certain boundary value problems -the so-called 
bounded quantum mechanical problems. 

Kuhn [13] presents considerable analytic and numerical detail 
obtained from expansions of the confluent hypergeometric functions 
in series of Bessel functions. His expansions can be obtained from the 
integral equation (15) by successive substitutions without any subse- 


quent rearrangement of terms except for collecting the J, and J 
terms. 


F 
Using a variety of methods, other expansions of the confluent 
hypergeometric functions have been found by Abramowitz |1], Buch- 
holz [4], Henrici [9], Karlin [12], Pignedoli [17], and Tricomi |19] 

\symptotic expansions of the confluent hypergeometric functions 
are given or referred to by Chang, Chu and O’Brien [6]. Such ex- 
pansions have also been considered by Moecklin [16] and Buchholz [5]. 

We wish to express our appreciation of many interesting and helpful 
discussions with Dr. G. I. Latta. 
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The Ergodic Theorem for Banach Spaces 


with Convex-Compactness 


VALENTINA KROTKOV and ISRAEL HALPERIN, F.R.S.C 


1. INTRODUCTION 


ET 7 be a bounded, linear operator on a real or complex Banach 
space B and let 


l 
n+ 1 


T 
T is said to be ergodic on B in the sense of J. von Neumann if for each 
fin B, T,,f converges to a limit when m becomes infinite 
In order that 7 be ergodic on B, the following conditions are neces 
sary: 


"f »() when n — @, for every f in B, 


AK tor all n, for some A 


Indeed, (1.1) follows from the relation 
(1.3) (1/n) T’ 


and (1.2) from the theorem that a convergent sequence of bounded, 
linear operators must be uniformly bounded [1, p. 80, théoréme 5] 
Clearly, if 7 is ergodic on B and 77,f converges to f*, say, then (1.1 


and the identity 
Pulit) T(7,f) 


together with the continuity of 7, show that 7f* ‘hia 

Because of (1.3), (1.2) implies that l/n)T"  < 3K tor all » and 
hence that (1/n)7"f will converge to O for all f if (1/n)7"f converges 
to zero for a set of f which span B, that is, whose linear combinations 
are dense in B 

Both (1.1) and (1.2) are implied by the condition 


(1.4) 7” | < K for all n, for some K < @, 


17 
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But if 7 is ergodic on B and (1.4) holds, then the identity 
ota ye) fe T"(T »f — f*) 


shows that 7),(7’f) converges to f*, irrespective of how the non- 
negative integer r varies, when n becomes infinite. 

In this note we discuss a method of Fk. Riesz [6; 7] to show that the 
conditions (1.1) and (1.2) are sufficient to imply that 7’ is ergodic on B 
if B is convex-compact. By definition, we shall say that a subset A of 
B (possibly all of B) is convex-compact if every decreasing sequence 
of non-empty, closed, bounded, convex subsets of A has a non-empty 


intersection. 


2. THe Riesz METHOD OF AVERAGES 


Suppose that 8 is an arbitrary Banach space and that 7° satisfies 
(1.1) and (1.2). Let S(f) denote the set of all ‘‘averages’’ of the form 


g=oftalf+...+67' 


where r is an arbitrary non-negative integer and the c; are non-negative 
scalars with sum equal to |. Riesz observed that for such g and for 


n ns 
- Tg ) ; 5 oes Pe es ee Tet apy. 


Hence (1.1) implies 


and therefore 


lim |7nf — Tnf| < lim |T7,/ T,g| + lim |Tyg - 


nN, Mex n n 


lim |g — Tmg| + lim [Tag — Taf 


m m-» 


2 lim |7ng — gi. 


Thus, following Riesz, 7 will be ergodic on B if there are elements g, 
in S(f) such that 


(2.1) lim |7ygi— § > 0 when 1 


Now there will be such g,; in S(f) if there is an element / in S(f), the 
closure of S(f), with 7h h. Indeed, there will then be g; in S(f) such 
that | g; > when 7 — ©; since 

DV n¥i ae Tale h)+(T,h —h) +h £1), 


the identity 7),h h together with (1.2) shows that (2.1) will hold. 
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Thus 7 will be ergodic on B if for every f, the mapping of S(f) into 
itself by 7 has a fix-point. 

Now let Sal be the set of elements wu in S(f) satisfying both 

u| <A\f|and| Tu — u l l)u | < 1/n. Then every S,(f 
is non-empty; for | 7,f| < K for all m by (1.2), and the identity 


T(Tnf) — Tnf | 7 <—| Ty] +s 


m-+ | m+] m+)" | 
together with (1.1) show that, for fixed n, 7),,f is contained in S,(/) for 
all sufficiently large m. It is easily verified that for fixed f, the S,(f 
are a non-ine reasing sequence ol bounded, ( losed, convex sets. If B is 


convex-compact then some element / will be in all S,(f); clearly, this 


h will be in S(f) and for every n, | Th h! <1/n so that Th h 


This proves: if B is convex-cempact then a bounded, linear operator 1 
1s ergodic on B uf and only if 7 satisfies (1.1) and (1.2 

Since metrically closed convex subsets of B are necessarily weakly 
closed (this is a theorem of Ascoli-Mazur, see Halperin, these ‘Trans 
actions, vol. 47, pp. 1-6 weak convex compactness Is equivalent to 
the strong convex-compactness as defined above. Hence our ergodic 
theorem includes the case of locally weakly compact Banach spaces 1S] 
and the equivalent case of reflexive Banach spaces [5], equivalent 
by the theorem of Eberlein. These writers assume (1.4) but their 
arguments are actually valid assuming only (1.1) and (1.2). The 
reflexive spaces are now known to include the uniformly convex 
Banach spaces of Clarkson [2], and in party ular, the L? spaces, Pp | 
We refer to [8] for references to other writers on this subject but 
parti ular reference should be made to the prool of an ergodic theorem 
7] 


by Riesz |7 


3. THe UNIFORMLY CONVEX BANACH SPACI 

We note that it is particularly easy to show that every uniformly 
convex B is convex-compact. Indeed, if S, are a decreasing sequence 
of non-empty, bounded, closed, convex sets, let y, infimum | g | for 
all g in S,, and let g, be selected from S, with un < | gn!) < un + 1/n 
Then the uw, are non-decreasing and bounded above since the S, are 
bounded. It follows that | g, | converges to a limit as n becomes 
infinite and the convexity of the S, together with the uniform con 
vexity of B implies that the g, converge to a limit h. This h is neces 


sarily in all S,, assumed closed 


1. ADDITIONAL (COMMENTS 


Call an element f in the Banach space T-ergodic if 7, f converges 


when » becomes infinite. Then (1.2) implies that the 7-ergodi 
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elements form a closed linear subspace of B. Our argument actually 
shows that f is 7-ergodic if S(f) is convex-compact, and that 7° is 
ergodic on B if S(f) is convex-compact for a set of f which span B, 
assuming (1.1) and (1.2). 

Finally, consider the special case that B is L': the space of f(P) 
summable on a set X of points P with a non-negative, countably 
additive measure | e| , and that 7’ is a bounded, linear operator of the 


type 


T{(P) = f(o(P)) 


for some mapping ¢(P) of X into itself. It is easy to see that L! is not 
convex-compact, except when X consists essentially of a finite number 
ot points. However, if X has finite measure and T is of this special 
type then for every measurable subset e the function 1,, which is 
defined to have value | if P is in e, and value 0 otherwise, has S(1,) 
convex-compact, in fact every sequence of elements g, from S(1,) 
has a weakly convergent subsequence. This follows from the fact that 
the g, are uniformly bounded (| g,(P) | <1 for all P) and X is of 
finite measure. Thus (1.1) and (1.2) imply that such a T is ergodic 
on such an L'(see [4]); but (1.1) may be omitted since it is auto- 
matically valid for all 1, when X is of finite measure, and the 1, span 
L': also, (1.2) can be deduced if | 7,1.!| < K!e! is known to be 
valid for all e, as pointed out by Dunford and Miller [3]. Thus the 
necessary and sufficient condition that a 7 of this type be ergodic on 
such an L' is that for some K < @ and for all m and all e, (see [3}]), 


‘(e)| +...+ | "(e)|) < Klel. 
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A Problem in Combinatorial Analysis 
N. S. MENDELSOHN 


Presented by R. L. JEFFERY, F.R.S.¢ 


INTRODUCTION 


gigantic in combinatorial analysis are usually considered 
solved when one can find either a difference equation or a recur- 
rence relationship which enables one to compute f(m) step by step. 
If it makes computation easier, an explicit formula for f(m) is a more 
desirable tvpe of solution. In common with other branches of number 
theory, certain published “‘solutions’’ of combinatorial problems 
merely replace the original problem by another one which may in 
fact be more difficult to solve than the original problem. The author 
of this paper has discovered in [1] and [2] a symbolic method of setting 
up and solving the recurrence relations for a wide class of combina- 
torial problems. The problem studied in this paper can be handled by 
the author's general method, but the computations are very laborious 
However, the formula, once obtained, can be verified much more 
easily and in this paper we will present only this easy verification. 
Numerical examples will illustrate the computational efficiency of 
our formula 


2. THe PROBLEM 


We will consider the following problem. Suppose we have a set of 
at+b+er4+... objects of which a are alike of one kind, 6 alike of 
a second kind, etc. Let C,la, b, c, ...] denote the number ofcombina 
tions taken r at a time of these objects. We will write C,[a’, 6, c, 


as an abbreviation for C,la, a, a, 6, b,c, ... |. Anexplicit formula for 


Ka 


“oo k 
C, fay GR viscnie) 


will be obtained. For the corresponding permutation problem, i.e., the 
problem of finding the number of permutations r at a time of the given 


class of objects, a recurrence formula will be obtained. Corresponding 


to the notation C,la, b, c, | we will use P,fa, d, c, | in the case 


of the permutation problem. 
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3. THe FUNDAMENTAL RECURRENCI 


‘first show that 


Cds Os) 62a tl 
satisfies the following recurrence relationship: 


mir 


(1 ) C,{a‘ ie a Kreg ae 2. ‘* las’. a’, ate ae ‘ 


lo prove this formula we isolate one set S of a, identical objects and 


observe how these objects enter into the combinations. The combina- 


tions which are counted by the expression 


ge} 
can be divided into the following mutually exclusive classes 
(() 


Those which contain no object of the set S. There are 


of these 


(1) Those which contain one object of the set S. There are 
of these 


(u) Those which contain uv objects of the set S. There are 


of these 


On adding the expressions obtained in (0), (1) 


wbtain formula (1) 


Formula (1) is quite efficient for computing purposes but an explicit 
expression will be given in the 
better 


next section which is actually much 


Pure Expriicit FORMULA 
\n explicit formula for 
‘ h Se ‘ 
( lay _ oe, segs 


can now be obtained. However, we will write down the formula for 


the case C,la', b/| only. A simple examination of this formula will 
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show how it may be extended to the general case. The formula in 


question is 


t+ | “ 

t+ l)m 
where the summation is taken over all integral values of the indices 
k, m satisfying the inequalities: 0 << k <1,0 <(m <j,0 < (a+ Dk 
+ (b+1)m <r. In practice, it will be seen that these inequalities 
contine k and m to a very small range. In the next section we will 
show that formula (I1) is quite easy to handle numerically. To prove 


formula (11), we first note that if 7 = 0, (IL) is replaced by the formula 


(111) 


summed over all values of & in the interval 


F 
QO<k < min (;, ) 
a - h. 


where LX] stands for the greatest integer less than ol equal to X 
When 7 1, <I11) becomes C,fa] i, if r <a, and C,[a}] = 0 fot 
r >a. This result is obviously true, so that (III) is true if 2 | 


Phe recurrence (1) for the case ¢ 1, Ry 1 is 
(IV) C,la‘| 


Assuming (II1) to be true when 7 is replaced by 7 1, and using 
(IV) one obtains that (II1) ts true for @ if it is true for 7 |. Hence 


(111) is true for all integral 2, by induction. Now C,la‘| C,la,;, b°| so 


that (11) is true when 7 0. By a similar application of (1) to the case 
t = 2 we can establish (II) by induction, to be true for all 7. This 
completes the proolt. 

One can obtain simpler formulae in case any of the a, are greater 
than or equal to r. In this case the number of combinations is actually 
independent of the exact value of a;. We will replace a; by the symbol 
© incase a; > rso that, for example, Cs5/2°, 3*, 54, 6°, 7°] will be written 


as C,{2°, 3?, '). The formula for C,la', 6’, 1 can be verified to be 


(V) C,la‘, b’, 07] = D> do (-1 
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where the summation indices k, m range over all integers satisfying: 
QO<ck<1,0Km<ej,0 <€ (atlket+(b+1)m <r. 
In the next section we will illustrate the efficacy of these formulae 


by numerical examples. 


5. SOME ILLUSTRATIVE EXAMPLES 
\s a first example we will consider the number of essentially different 
l1-card canasta hands which can be formed from a 104-card pack. 
(We do not count the red threes as part of the pack since these are 
always replaced in a hand.) The pack consists of four black threes, 
four jokers, eight aces, eight twos, eight fours, ...., eight kings. 
Cards of a given denomination are considered equivalent regardless 
of suit so that the number of ways of choosing eleven cards from such 

a pack is C,,[8'’, 42]. Substituting in (II) we obtain 


= us 12\/ 2\/24 — 9k — 5m 
‘ gl 9 s rem( 4 bi & ) 
( ri ’ i | > p> \ l) k m 1l — 9k — 5m 


summed over all pairs of integers (Rk, m) subject to the inequalities 
0O<¢g¢k 12,0 Cm < 2, 0 < 9k + 5m < 11. The pairs (k, m) satis- 
fying these inequalities are: (0,0), (0, 1), (0, 2), (1, 0). Hence, 


: » ( . 
Ci:[8", 4°] = i“ 2( 1°) + 14 (1) = 2,440,634 


It will be noted that our solution requires the computation of four 


terms only. A solution based on the partitions of eleven would require 
the computation of 52 terms each of which is a product of binomial 
coefficients. A solution by J. P. Ballantyne [8] would require the con- 
struction of a table containing 12 rows and 16 columns, 1.e., one 
containing 192 entries. 

As a second example, we consider the problem of finding the number 
of essentially different bridge hands. By this we mean the number of 
different bridge hands which can be formed when the cards in each 
suit from 2 to 9 inclusive are ranked equally. The required number is 
Ci3/84, 17°). By (ID) 


fe age - 4 bag be Ok - 
seashell X »» : m 13 — 9k — 2m 


summed over all pairs (k, m) such that OC k <4, 9 Om < 20, 
0 < 9k + 2m < 13. Hence, the admissible (k, m) are: (0, 0), (0,1), 
(0, 3), (0, 4), (0, 5), (0, 6), (1, 0), (1, 1), (1, 2). This yields 
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Cis(8', 12°] (=) - 20( 3!) n i90(3”) 1i40( 2”) is1s(°8) 
iss0i( 2°) + 38760( >!) "id n so(?°) 760 


34,033,880. 


Here again, the saving in computation is very considerable when 
compared with the usual counting methods. This result is in agreement 
] 


with that given by Ballantyne in [3] and affords a comparison between 


the two methods 
6. THE CORRESPONDING PERMUTATION PROBLEM 


We will denote by P,la, ), c, | the number of permutations taken 
ratatime of a+6b+cct. objects of which a are alike of one 


kind, 6 are alike of another kind, ete. We will use the notation 


and P,fa’, 7] 


to represent the permutation analogues of 


C,[ai', a2*,...,a4'] and C,[a‘, 7]. 


We will confine ourselves mainly to the problem of finding the permu 


tation analogue of (1) namely 


min(7,a@e 
(XT a r A r Ke 
(VI) P,lat,ai,...,0¥%)= > (") p, Leb as. 
u 


0 
‘ 


Before proving (VI) we note that the formulae 


PAV r)! 


are easily obtained from (VI) by induction. These formulae are, in 
fact, classical results 

The proof of (VI) is analogous to the proof of (1). We isolate one 
set S of a, identical objects and then divide the totality of permu 
tations into the following exclusive classes: kor v = 0, 1, 2, 3, 
min(r, a;), we form the class P,,, of all permutations containing 
exactly v elements of the set S. Then, obviously 

min(r,a¢) 


2 ks k 
PGs (Oe pens phe | bs P,, 


To compute P;,, we note that each permutation in the set P,,, can 


be obtained by forming a permutation from r — v objects chosen from 
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the totality of objects except those in S and then introducing into 
this permutation v objects taken from S. The number of ways of 


introducing these v objects into the permutation is 


(’), 


and the number of permutations into which these objects must be 


introduced is 


Hen Cc, 
P 


This completes the proof of formula (V1 
\s yet, the author has been unable to transform the recurrence (VI) 
gives rise to a 


In spite of this, formula (V1 


into an explicit formula 
rapid computational procedure, and in spec ial cases does lead to 


expli¢ it formulae 
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Note on a Combinatorial Formula of Mendelsohn 


LEO MOSER 


Presented by R. L. JEFFERY, F.R.S.C. 


HE object of this note is to give a short, direct proof of result (11 


of the preceding paper. The result in question is 


(II) C,[a'b4} = YY (- NG) 


(’ titr—1l1—(@ti1)k—- (04 = 

r (a+ 1)k {b+ 1) m 
where the summation ranges over all values of & and m for which the 
binomial coefficients involved are positive 


Using the set up of Ballantyne’s paper [3], we have C,[a', b/] is the 


coefficient of x" in 


f(x) 


Hence 


Cra’) = ED (-yee(? 


where ¢ is determined by r (a + 1)k 


This completes the prool ol (11). 
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Some Remarks on Laplace’s Method 


P. G. ROONEY 


Presented by M. WYMAN, F.R.S.¢ 


| i his book, The Laplace Transform, D. V. Widder finds asymptotic 
formulae for integrals of the form 


| exp (kh(x)) o(x) dx 
ea 


as k—» , under the conditions that A(x) has a simple flat maximum 
at a point 4 c between a and A, i.e. h’(c) 0, A’ (c) < 0, and (x 
has, at worst, a simple zero at x ¢ and does not change sign there, 
together with certain other conditions 

These results can easily be generalized to the case in which h(x) has 
a higher order maximum at x c, and @(x) has a zero of arbitrary 
order there, though (x) is still not allowed to change sign at 4 ( 
The generalization is contained in the following two theorems 

The first theorem deals with the case in which the maximum of 
h(x) occurs at one end point of the interval, and the second theorem 
with the case in which the maximum of A(x) occurs at an interior point 
of the interval 


THEOREM | 


, h(x) is non-increasing in a 
a)’ o(x) L(a, bh), o(at+) exists and o(a4 


then 


‘a 
| exp (RA(x)) (x — a) (x) dx 
eva 


" v+l a : 
~ (—(2m)!/kh (a)) (p(a + )/2m) exp (kh(a))P((v + 1 
as k 
Proof. We may assume ¢(a+) > 0. Then for every ¢ such that 
0<€« < d(at+) andO < e€ « h°*"’ (a), 6 exists such that 


29 
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bla ) < | =e ola+) + €, 


IR(h (x 


i, 


l 


ah 
exp [Rk(h(a + 4) h(a))} | o(x)\ (x a) dx. 
wa 


Since h(a + 6) h(a) < 0, it is clear that J,’ = O(a‘), 0 
Using (i) and Taylor’s theorem with remainder, we have 


A, | exp [RA (E) (x ay” /(2m)'!}(x — a)’ dx < I,’ 
e 


ea 
As | exp [RA (E) (x a)" /(2m)'} (x — a)’ dx, 
va 


where a < & <a + 6. From (ii) we only strengthen these inequalities 
by replacing h@” (&) by h@” (a) e on the right and by hh” (a) + « 
on the left. Thus 


+6 


ea 
A, | exp [RBy(x a)” /(2m)' (x 
e a 


A» | exp |RBo(x 


In the left hand side of these inequalities let 

u RBy(x a)” /(2m)!, 
and in the right hand side, let 

u kBs(x — a)""/(2m)!. 


Thus the inequalities become 
Ke 


. » 
, v+1) /2m ; +1) /2m)—1 
(—(2m)!/RB,) (A, /2m) | du < I,’ 
Jo 
ek, 
v+l 2 u v+l tn 1 
(2m)'/RBo) A»/2m) | e ou du 
et 


where C, B,6"/(Qm)! (1 ey 
Thus, 


v+l 2 . - * 
(2m)!/B,) °(A,/2m) T((v + 1)/2m) < lim 7,’ 


v+l 2m ‘ . i ‘ 
(2m)!/Bo) (Ao/2m) T'((v + 1)/2m). 
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Hence, since € is arbitrary, 


. v+l . v+l 
lim k lim R 

( "(bla + 
COROLLARY | 


(he <8 
2) A(x C*"(b— 9 

0, h(x) ts non-decreasing 
Lia, b), d(b 


ina 


O(x) exists and (b 


then, 
/ 
| exp (Rh(x))(d v) (x) da 
e a 
2m) exp (kh(b)) P((v + 1)/2m) 


™~ | (2m)!/kh su |: (b(b 


COROLLARY 2. /f (1), (2), and (3) of Theorem 1 hold, with the excep- 
0), then 


tion that d(a 


ef 
exp (Rh(x)) (a a)” d(x) dx o(k *""™ exp (Rh(a)) 
va 


THEOREM 2. // 

(l)ha <b n<b<b+yn ; 
(2) h(x) C2" (b nqx« t+ n), A (b ee 
2m 1), AY™ (bd) Q), h(x) is non-decreasing in a s < ), non 


increasing in b as, 


(3) ix bh’ o(x) Lia, c), 


(b+), ob exist and are not 


then 


| exp (Rh(x))\x — b) (x) dx 
va 


~(—(2m)!/Rh°™ (b)) )/2m) 


v+l 2m { 
((o(b-4 ) + o(b 


exp (kh(b)) T'((v + 1 


as k 
Proof. li we write the integral as the integral from a to 6 plus the 
integral from 6 to c, and apply Corollary | to the first and Fheorem | 


to the second integral, we arrive at the result directly. 


As an example, we may find the asymptotic value of 


el 
kr* 6 
| e sin x dx 
e 0 
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This integral fulfills all the hypothesis of Theorem 1, since sin® x may 
be considered as x* @(x) with (0+) 1. Thus we have 


al 
kr* a 7 ssa 
| e sin xdx~ gk I (5) 
J 0 


lurther, since 
er . 
z* - 6 k~1 . 
| e”” sin xdx| <e | ¢6 ds, & > I, 
ea) e/ | 


which is of smaller order than k-7/* as k — ©, we can say that 


| e** sin’ x dx ~ ik” ri) as k 


Widder's results may be generalized in a different direction. The 
integral may be a double integral, and 4 and ¢ functions of two variables. 


he generalization is contained in the following theorem. 


(1) Ris a measurable region of the x — y plane, such that 1/ 
Ss Px: apr ly Bb) < 6}, 
then for some 69, Ss, CR, 
(2) d(x, y) C*(Sa.) (\ L(R): 6. = @, Os ‘Gas? bro, < O, 
bry < 0, at (a, b); for each 6,0 <6 < by, O(x, vy) < M(6) < g(a, b) for 
(x, ¥) R Ss, 


(3) W(x, y) L(R), w(x, y) is continuous at (a,b) and Pla, b) ¥ 0, 
then 


| | exp (k(x, y)) W(x, y) dxdy 
. JR 


~ 2m la, b) exp (k d(a, b))/R( dr: (a, 0) by, (a, 6 dry (a, b)) 
as k— @, 


Proof. We may assume (a, b) > 0. Let 


a Ecos 0 nsin@, y Esin 6 + n cos 6 
where 
tan 20 = 2¢,,(a, b)/(¢,,(a, b) — 9,,(a, b)); 
let 
d(x, vy) P(t, n), w(x, y) W(E, ). 
Then since 
O(x, Vv) 


O(€g, n) 
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xp tR( d(x, vy) — ofa, b))} W(x, y)dxdy 


ee 
“4 


f exp (k((E, 7) — &(0,0)) WE, 1) d&dn. 
Os 


By 2 and since ;;(0,0) <0, ,,(0,0) <0, :,(0,0) = 0, we 
can, for every «, 0 <¢€< min. (—%;, —®,,, V(0,0)), choose an 


ale), O <a « (460)?, such that, 
(i) Ay P::(0, O) e< Per(t,n) < By ?::(0,0) + € <0, 
e< ;,( +7 
A, = ®,,(0, 0) < ,,(&, ) < o = $,,(0,0) + « < 0, 
(i) 0O< Cy Vv (0, 0) < gE, 6. WV(0,0) + «, 
for all (&, 7), a<cé< 
Let R, { (x, y) ax<t<a; <<a}. Then 


. . | 
ai iJ Jr. J iF R,§ xP tk(o(x, y) — O(a, 6))} 


Consider first /,’’. Since, for (x, y) R R., (a a)? + 


ve have, by (2), 


0 « | ig S exp a Ma) Pla, )) if s V(x, y) dxdy 


oo. that ,’’ w*) where 0 < w < 1, since o(a,b) > Mila), and 


v(x, y) 
Now 


yy) g(a, b))} wir, y) dxdy 


1) b(0,0))} W(E, 1) dédn 
" [" ’ [ fg | 
+ + 
. ) et e ¥ e ae m 
exp 1k( P(E, n (0, 0))} WE, n) didn 
+ J. 
By (2) and Taylor’s theorem, 
(0,0) = 4( ae. (E, On) E + 24 ¢, (OE, On) En + ,, (OE, On) n ), 
one <3 


Consider J;. By (1) and (it), 
C\ | | exp (4k(AiE Jetn + Aon )) dtdn < J, 


C» | | exp (4R( Bie +t Jetn + Bon’ )) dédn. 
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In the left hand side, let u SRA \&?, v? = kA’, and in the 
right hand side, let u? SRB, &, v? ‘kB.»?®. Then the in- 


equalities become 


oh D eh dD, 
(2C,/k(A,A>2)') | | exp ( (u” + 2€(A,Ao)* uv + 0°)) dudv 


eh Db eh De 
(2C./k(B,B2)’) | 


exp(— (u° 2¢(B,B,) ‘uv + v°)) dudv, 


where D, (—4A,)?, (1 1,2) and D, a( 1B,)), (1 
Hence 


I 


(2C,/(A1As)’") | | exp ( (u” + 2€(A,Ao) *uv + v°)) dudv 
eit ei 


lim kJ; < (2C2/(B,B:)!) | 
exp ( - (n° + 2e(B,B.)* uv + 0°)) dudv 


and thus, since ¢ is arbitrary, and the integrals in question converge 
. i 
uniformly in ¢ for 0 < € < e«) < $(B,B.)?, we have 


ium kJ, (2 (0, 0) / (de. (0, 0) &,, (0, 0))! | | eo” *"” dudv 
kw» et eit) 


w Wa, b)/2(¢,,(a, b) o,,(a, dD) br, (a 
Similarly, 
lim kJ» lim kJ; lim kJ, nr (a, b)/2(¢,,(a, b) $,,(a, b) 


A 


bry (a 
Hence 


lim RJ, lim k/,’ lim kR( J; + Jo + J3 + J) 


2m (a, b)/(drr(a, b) b,,(a, 6) — br, (a 
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A Remark on Curves of Order n in n-space 


P. SCHERK, F.R.S.C 


1. Let R, denote real projec tive m-space. An arc A [a curve C] is 
the continuous image of an interval fof a cire le] in R,. The order of A 
is the upper limit of the number of points which A may have in 
common with any [linear] ( 1)-space. We exclude the case that A 
degenerates into a single point. Then the order of A is not less than n 
Let A” {c"| denote an arc {a curve] of order 1 in R,,. It necessarily is a 
Jordan are {a Jordan curve]. 

Let s,s, denote points on the are A. We call the point » differ 
entiable if it possesses osculating subspaces L,"(s) of every dimension 
[p 1, O, I,..., mf. Let Loi"(s) be empty space. Suppose we 
have defined L,"(s) and postulated its existence. Then we require 
that the (p + 1)-spaces through L,"(s) and a point s’ converging to s 
have a unique limit space. We call it the osculating (p + 1)-space 
Loyi"(s). Thus Lo"(s) is the point s itself and L,"(s) is the whole R, 
We call A differentiable if each of its points is 

It is well known [1] that the end-points of an arc A” are diffet 
entiable. Thus, the interior points of A" are “‘one-sidedly differentiable.”’ 
This implies that A" is rectifiable and that it is differentiable every 
where with not more than a countable number of points as exceptions 

2. In a recent paper [4], Schoenberg proved with analytic methods 
that every C" can be uniformly approximated by analytic C"'s if ” is 
even. He also proved for odd n’s that every A" can be uniformly 
approximated by ares of analytic C"’s with the same end-points. In 
this note, a geometrical proof of the following weaker theorem will be 
outlined 

Every C" can be uniformly approximated by differentiable C"’ 

Here m may be any integer greater than one. Our proof is based on 
(2] and [3]. The author should like to acknowledge the stimulus he 
received from Dr. Schoenberg and Dr. Motzkin 

3. Given an (n — 1)-space / and an interior point s of an are A’ 
Let A” be a small closed subare of A’ which contains s in its interior 


and which meets / nowhere else. We then call s a point of support 
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(of intersection] with respect to £ if the end-points of A” lie on the 


same side of E [are separated by /]. The following lemma is quoted 
from [3]: Suppose s is differentiable and - contains L,"(s) but not 
Lpyi"(s). If p is odd [even], then s is a point of support [of intersection] 
with respect to EF. 


, 


1. Let s and s’ be two differentiable points on a curve C"; s # 5’. 


‘n 


They divide C” into two arcs A" and A’". 


n 


The osculating spaces L,"(s) and Ly ~p—2"(s’) are readily seen to 
span an (” — 1)-space which does not meet C” outside of s and s’ and 
which does not contain L,4i)"(s) or Ly_-pia"(s’) [-1 < p cn — II. 
Altogether there are n + 1 subspaces of this kind. We verify by in- 
duction that they subdivide R, into 2” simplexes. From the above, 
A" {A’"| lies entirely in one of these simplexes. The lemma quoted in 
$3 implies that the simplex S that A’” lies in is determined by A” 
Let A’ be a third are of order m with the same end-points s and 5’. 
Suppose A" U A’’" is differentiable at s and at s’ and shares with C” 


/ 
2"(s°) 


the property that the (nm — 1)-space through L,"(s) and L,_, 
supports or intersects A" \U A’’" at s if p is even respectively odd 
[p= 0,1,2,...,% 1]. Then, from [2] and [3], A" U A’’” is a curve 
of order n. Krom the above, A’’" will also lie in S. There are differ- 
entiable arcs A’’” with the required properties. As a matter of fact, 
A’’"™ can be an arc of a rational curve of degree n. 

5. Given the curve C" of order n. We introduce an elliptic metric 


"" 


into R, and decompose C" by means of differentiable points into 
sufficiently small subarcs A’”". If our decomposition is fine enough, 
then each of the simplexes S determined by the A’” will be arbitrarily 
small. Replacing consecutively each A’ by a differentiable A’’", we 
obtain a curve of order n which is differentiable everywhere and which 
lies in the union of the simplexes S. This proves our theorem. 
Schoenberg's theorem could be proved in this fashion if the fol 
lowing conjecture is true: To any finite set of differentiable points on a 
C" there exists an analytic curve of order n through them which has 


at each of them the same osculating spaces that C” has. 
REFERENCES 
J. Hye_MsLev, Introduction a la théorie des suites monotones, D. Kgl. Danske 
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The Application of Fourier Transforms 
in Physical Problems 


W. H. WATSON, F.R.S.C 


HE application of Fourier transforms, and integral transforms in 
pda, is of quite modern interest for it is still evolving. Its 
study exposes in a fundamental way the interactions of mathematics 
and physics. In the main, books about it suffer from one or other of 
two faults. Mathematical interest is usually confined to the theory of 
some special transform such as Mellin or Laplace, and the applications 
discussed are often trivial, illustrating basic points of theory it is true, 
but nevertheless falling short as a guide to the detailed analytical 
technique needed to handle real problems. This afternoon it would be 
out of place to venture on mathematical elaboration that properly 
belongs to sustained study of the subject. It is appropriate, however, 
to recall in mathematical form the basis of integral transforms. 

Given a function A (a, x) of two variables a and x, and the function 
f(x), the function F(a) defined by 


F(a) = | f(x) K(a, x) dx 
e 0 


is, When the integral is convergent, the integral transform of the func- 
tion f(x) by the kernel K (a, x) 
In the Laplace transform K 
In the Mellin transform AK 
In the Fourier cosine transform K (a, x COS ax 
In the Fourier sine transform A sin ax 
In the Hankel transform of order v A ,x) = (ax)? J,(ax) 
In each case there is an appropriate inversion formula by which 
f(x) may be deduced from F(a) {7}. 


The Fourier transform of f(x) is 


(1) f(a) | f(x) e* dx 


and the inversion relation is 
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) 


Ll (*; a 
(2) f(x) f(a) e * da 
Te , 


f(x} must satisfy Dirichlet’s conditions and 


| f(x)| dx 


It is forty years since Sommerfeld applied Fourier transformation to 


must converge. 


treat the propagation of a signal in a dispersive medium. The problem 
arose from the objection offered to [Einstein's Special Relativity theory 


an objection supported by citing examples of media in which, since 


dV/dX 1s negative, the group velocity exceeds c, and put forward in 


the belief that the group velocity is necessarily equivalent to the velocity 
of energy propagation. The particular problem considered by Sommer 
feld [8] is that of a signal which originates at z 0, at instant ¢ 0, 
and continues indefinitely thereafter as a harmonic oscillation of 
frequency wo (radians 

Phat the Fourier integral of a non-terminating harmonic wave does 
not converge may be seen from equation (1). The difficulty may be 
circumvented by deforming the path of integration into the complex 
domain of the variable wy). For example, for 


ae 


e 


| 
) 


“211 (wy wW) 


3) 


provided Im (w 


ot 
(4) ; m ; dt. 
] 27 (Wyo Ww) 

This integral is evaluated by the theory of residues, the contour being 
closed by a semicircle at infinity on which, by Jordan’s lemma, the 
contribution to the contour integral is zero. When t > 0, the path ol 
integration lies in the negative half-plane at infinity thus enclosing the 
pole w wy. Whereas when ¢ < 0 it lies above the real axis and excludes 
the pole, as shown in Figure 1. This gives the essential discontinuity 
in the form f/(t). 


Po obtain the Fourier transform of 


we eV idently have to combine 
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l r ‘ 
at aa dt and 
the 


> 0 tor the first and Im(w 1) < for 


requiring Im(w — w 
second. It is in this sense that we write 


and regard 6(u as the transform of 


| 
. é 
25 


6 has been called the Dirae 6-function. Schwartz 16] calls it the Dirac 


measure. It is not a function. | shall return to this topic later in the 


lecture 


Phe broad use of integral transforms to solve differential equations 


has developed from Heaviside’s operator calculus, devised by him to 


solve transient problems in physics and electrical cngineering The 


Important core of Heaviside’s method is that it corresponds to the 


possibility of computation with numbers instead of analysis because 
it starts with what is given, and evolves the answer explicitly by direct 


operations on the data 
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If one wished to characterize the role of Fourier transformation in 
the theory of physical systems, it may, I think, be fairly stated that 
fourier transformation is the natural extension of Fourier series to 
unbounded spaces permitting the representation of discontinuous 
distributions in space and time by methods conformable to analysis. 
Whether discontinuity is a physical necessity does not really matter: 
it is mathematically simpler than the representation of the actual 
rapid transition according to any special law that might be physically 
appropriate, but which, nevertheless, is judged irrelevant to the main 
idealized purpose of the calculation. This use of discontinuity in 
representing nature is quite analogous to the use of impulses in dynam- 
ics. 

It has already been pointed out how the switching on of an oscillation 
is managed by treating Fourier transformations in the complex domain 
of the transform variable: a great deal of the subtlety of calculation 
emerges from this analytical way of representing discontinuous distri- 


butions. Since the physicist expects to see in the mathematical develop- 


ment forms of which he recognizes the physical relevance, he will tend 
to by-pass whatever he thinks he may dismiss as mathematical ritual. 
This is in some respects very unfortunate, because he is often thereby 
deprived of the stimulus to appreciate the delicacy of the conceptions 
being treated and so may not learn to use the very sharp tools of the 
mathematician. 

Let us now make a list of topics in the study of which the application 
of Fourier transformation has piayed an important role for theoretical 
physics. 

(a) In the theory of vibrations, Fourier series are the main tool, 
integral transforms apply usually only to the time dependence. 

(b) The diffusion equation and that of heat conduction we might 
regard as simplified versions of the transport equation and indeed as 
representing the whole field of continuous linear stochastic theory. 

(c) In elasticity stress systems due to discontinuous pressure dis- 
tributions have been treated 

(d) Hydrodynamical problems involving discontinuous boundary 
conditions. 

(e) Electromagnetic theory and electromechanics have been the 
field of many successful applications. 

(1) kourier transformation is peculiarly suited for the study of wave 
equations in general and especially to provide powerful methods for 
handling singularities and discontinuities. 

(g) In modern physics one can point to the use of wave functions in 


momentum space to treat scattering problems, the use of Mellin and 
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Laplace transforms in the cascade theories for cosmic rays, and the 
wide use of Fourier transformation in quantum electrodynamics and 
wave field theory. 

For the mathematician, however, the application of Fourier trans- 
formation to physical problems consists of its use to solve (i) partial 
differential equations and systems of them, (ii) integral equations and 
systems thereof and in general, (iii) integro-differential equation 
systems. Of course, Fourier transformation has its role in mathematics 
itself, and it is an important tool in statistics. 

This is an impressive list and naturally not much of it can be exhibited 
in this lecture. | propose to pick out a few topics to illustrate what has 
been going on. 

The theory of the slowing-down and diffusion of neutrons has 
received much attention because of its importance for understanding 
nuclear reactors. We shall consider a very small part of the theory of 
neutron diffusion under extremely rigorous (indeed unrealistic) physical 
assumptions designed to simplify the problem [5]. We begin with the 
transport equation by which the rate of change of the neutron density 
following the motion is accounted for by capture, scattering and the 


presence olf neutron sources: 


dy(r, B, t) 
dt 


P(B)W(r, B, t) + f Paya, 8’) dp’ yp(r, B’, t) 


with 


d y Oy 


(6 
" dt ot 


+ div (vw). 
The functions ¥ and q denote respectively the density of neutrons and 
of neutron sources both per unit volume at the space position ft at 
time ¢, and per unit range of the symbol 6 that stands for all other 
parameters required to specify a neutron as, for instance, its velocity V 
or direction of motion S. The probabilities per unit time that a neutron 
is scattered or captured from the range d8 of 8 are respectively P,(8) 
or P.(8). Their sum is P(8), the total probability of collision. The 
relative probability that a neutron originally with parameters 8’ will 
be scattered into the range (8’, B’ + dB’) is denoted by f(’, B) dp’. 

The simplifications we shall introduce are: 

(a) scattering without change of energy, all neutrons having the 
same speed v which will be taken as unit; 

(b) scattering and sources ISsSOtrOpIc ; 


(c) stationary state 
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The mean free path (/) being taken as the unit of length, the unit of 
time is //v. Then P L, Pea, ane ?, 1 —a@ and we have in 
place ol (6) 

(7) S.grad y(r,S) + p(r, S) py Hole 
T 


where 
Wo (r J vir, $1) as’ 


is the neutron density at r, and dS is the element of solid angle 
The integro differential equation (7) to determine y is equivalent to 


an algebraic equation on Fourier transformation, viz. 
. ‘ . l a j i 
(3) ik. S o(k, S) t o(k.S) | oo (k) t ho (k) 
T T 


where 
o(k, S) y(r, S) exp (ck. r) r drsin@d0d@ 
. 


and (kK) and Ag(K) are the corresponding integrals for po(r) and 

go(r). The differential operator S.grad becomes the scalar product. 
Now rearrange equation (8) thus 

, ( 

bar (1 8) 


Integrate over all directions S to obtain 


o(k, S) | a) oy (k) t A (k i. 


(9) go(k) = {(1 — a) do(k) + Ao(k)} (,* Lr 


ik. S° 
Since the integral is equal to Rk! tan>! k, we may obtain the transform 
oy (kK) of the neutron density directly in terms of the transform A9(k 


of the density of sources: 


1 
t k 
(10) (k) = — do(k 
R (| a)tan ‘k 


which has the form 


(11 go(k) = g(k) do(k). 


Because the operations multiplication and convolution correspond 
on transformation and since g(k) is the Fourier transform ot e~', 4rr’, 


the equation (11) reads after inversion 


. k 
e , ’ 
(12 Yo(r Pyar’, 


i Gol 
inR- * 


where R 4 r’ , showing how the source element at r’ contributes 


to the neutron density at r. The symbol dr’ denotes the volume element 
in the space ol Tr. 
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As a second problem we shall consider the decay of resonance radia- 
tion by spontaneous emission to determine the shape of the resonance 
line. The physical system consists of one atom initially excited and the 
radiation field initially devoid of photons. It is represented by a state 
vector w which is a function of the state of the atom and the numbers 
N, of the photons in the states ¢. Since we are dealing with resonance 
radiation we restrict ourselves to the excited state A and the ground 


state B of the atom and introduce the two relevant basis vectors 
a(t) ¥(A, 0,) and A(y, ft) V(B, i). 


The fundamental quantum equation 


hoa y 


3 
(13) i OL 


Ty 
is expressed by the pair 
da 
at 
oh 
oe 
at 


| W(v) b(v, t) S(v) dv 
J 

* 
Vv at 2arilv v)b 


where AW (vy) is the interaction operator of the Hamiltonian v being the 
photon frequency, and 
Sav | 


( 


hv, Eu Es, S(v) 


and V is the volume of the hohiraum. Equation system (14) is to be 
solved subject to a(Q) 1, bv, 0) 0. On Fourier transformation 
with respect to ¢, the transform variable being z (and W* the complex 
conjugate of W) we obtain 


(15a) 2da(z) | W(v) b(v, 2) S(v) dv + a(0) 
J 0 


~ * 
(15b) sh(v, z) 2riv — 1) b(v, 2) + W a(z). 


Solve (15b) for b(v, 2), substitute in (a) and fiad a(z), which in turn is 
used to give b(v, z). On inversion we obtain an integral formula for 
b(v, t), a tormal solution from which the shape of the resonance line 
has been worked out [4]. 

What | have stated regarding the special connexion of Fouriet 
transformation with discontinuity applies in a peculiarly significant 
way to the treatment of waves. 

For Laplace's equation Al Q and for the iterated Laplace's 
equations A"/ 0, which are of elliptic type, there are no solutions 
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other than the analytic functions. However, for a hyperbolic equation 


such as 


al aU 
Ox P ay” 


= () 


the general solution will be 


/ A(x + y) + B(x — y) 


where A and B may be continuous functions, not necessarily differen- 
tiable, discontinuous functions and even not functions at all in the 
sense of classical analysis. The introduction of such discontinuous 
solutions in the theory of partial differential equations is of great 
importance to physics. [It is remarkable that substantially the un- 
orthodoxy initiated by Heaviside is still continuing. The latest adjust- 
ment of mathematics to physical intuition is exhibited in the note- 
worthy Theory of Distributions due to Schwartz. His book on this 
subject appeared in 1950 and members of this audience who attended 
the meeting of the Canadian Mathematical Congress at Vancouver in 
1949 will recall the excellent series of lectures he delivered on distribu- 
tion theory. This theory has found a proper mathematical basis for 
handling discontinuous solutions for partial differential equations and 
of escaping from the mathematical category of functions to treat 
Fourier transformation. 

lor the physicist a constant harmonic wave of infinite duration has 
an infinitely sharp line spectrum. It is represented by a 6-measure 
The inverse relation surely requires the Dirac measure 6 to have a 
Fourier transform. Here physical intuition has kept ahead of mathe- 
matics. For 
F 5(w wo) b(w we dw e 
. 


2n 2r 


Now this formal equation we have written down does not fall under 
the terms of classical analysis, in fact in Lebesgue’s theory of integration 
the integral is zero. Accordingly this formalisra can be kept ‘on the 
rails’ only by physicat understanding, whereas the value of mathemat 
ics is that logically correct answers are assured by obeving the rules of 
calculation, Schwartz's work seems likely to improve the situation but 
there will be a period when notation will be muddled —a matter of 
serious practical effect in modern theoretical physics. 

It is ten years since Heisenberg [3] invented the S-matrix that plays 
so important a part in quantum electrodynamics and the theory of 
wave fields. In his papers he represented the asymptotic wave functions, 


representing for example emergent particles from a reaction, by Fourier 
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transforms. For instance 6,(ky — Rk) associated with time dependence 


exp (—tkoct) represents an outgoing spherical wave as may be seen as 
. 
J 


| Lf ie ¢: ) 
(27) | (ky — k) ak. 


We assume r large: since ky is real Im k < 0, and hence the main 


follows. 


elr e Rrs . 
doddtk dk 
J 2m (Ro k) 1g di 


(Im(ko k) 


contribution to the integral is 


] 
(29r)° 
which in turn comes from the vicinity of ko, the result being 


koe 
lar 

In the quantum electrodynamics of Schwinger, Feynman, Dyson 
et al., one has to deal with the singular wave functions associated with 
sources and sinks to represent the creation and annihilation of particles 
Quite apart from the question of relativistic invariance, the essential 
mathematical difficulties are concerned with the fact that the elemen 
tary solutions of the various wave equations (due to point singularities) 
are properly speaking not functions at all. Nevertheless one writes 
down expressions for their Fourier transforms. 

Fourier transformation enters as the obvious means for dealing with 
the convolution integrals that appear in the calculation of the wave 
functions representing the quantum stochastic process [1]. For example, 
a particle starts from A travels to B where it is scattered, what is the 
probability amplitude to reach C from A after such a single scattering 


process? It has the torm 


, , 


dx 


J vx x’) f(x 
e 


and its Fourier transform is G(k)F(k) which is a. straightforward 
product. The transforms are all singular and it is in the handling 
of them by conventional means that some of the difficult divergencies 
in the theory arise. In a recent paper, Giittinger [2] has shown how 
Schwartz's Distribution Analysis automatically leaves a place for 
mathematical operations equivalent in effect to renormalization, 
regularization, or cut-off procedures: it does this quite naturally and 


brings to light that the theory as at present conceived is physically 
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incomplete. There is a good expectation that distribution analysis may 
remove some of the mathematical difficulties of quantum field theory 
especially as applied to mesons 


Thus we view today an interesting prospect. New mathematical 


tools have become available to the theoretical physic ist, but perhaps 


even more promising is the realization that mathematicians have been 
too ready in the past to impose on physics the unnecessarily restrictive 
categories of their own inventions. It is significant that the application 
of Fourier transforms has been the source of these developments that 


illuminate the role of mathematics in relation to physics. 
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